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Abstract 



We study the following ID two-species reaction diffusion model : there 
' is a small concentration of B-particles with diffusion constant Db in an ho- 

■ mogenous background of W-particles with diffusion constant Dw, two W- 

, particles of the majority species either coagulate {W + W — > W) or anni- 

^ ! hilate {W + W — > 0) with the respective probabilities Pc = — 2)/(q — 1) 

\ and Pa = — 1); a B-particle and a W-particle annihilate {W + B — > 0) 

with probability 1. The exponent 9{q,X = Ds/Dy/) describing the asymp- 
totic time decay of the minority B-species concentration can be viewed as 
a generalization of the exponent of persistent spins in the zero-temperature 
\ Glauber dynamics of the ID (/-state Potts model starting from a random ini- 

Q \ tial condition : the W-particles represent domain walls, and the exponent 

J-^ ' 9{q, A) characterizes the time decay of the probability that a diffusive "spec- 

tator" does not meet a domain wall up to time t. We extend the methods 
^ . introduced by Derrida, Hakim and Pasquier {Phys. Rev. Lett. 75 751 (1995); 

Saclay preprint T96/013, to appear in J. Stat. Phys. (1996)) for the problem 
of persistent spins, to compute the exponent 9{q,X) in perturbation at first 
order in (g — 1) for arbitrary A and at first order in A for arbitrary q. 

05.40.-hj, 02.50.-r, 82.20.-w 
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Typeset using REVT^ 



I. INTRODUCTION 



The one- dimensional Ising or Potts model evolving according to zero-temperature 
Glauber dynamics from a random initial condition is one of the simplest systems for 
which domain coarsening [0] can be studied in great details. The possibility of writing 
closed kinetic equations for the expectation value of each spin and for the equal-time two- 
point correlation functions can be used to obtain various exact results, such as the 
growth in time of the characteristic length of the coarsening like t^^^, as expected in general 
when the order parameter is not conserved. More recently, it was shown that more refined 
quantities such as the fraction of spins which have never flipped up to time t P,^, or the 
distribution of domain sizes could also be studied by mapping the problem to an exactly 
soluble one-species coagulation model {A + A — > A). 

The zero-temperature Glauber dynamics of the g-state Potts model starting from a 
random initial condition is related to various reaction-diffusion problems. The simplest 
relation deals with the dynamics of domain walls W ( |^ and references therein), that 
diffuse and react whenever they meet according to 

W + W — > W coagulation with probability pc = (|5f 



W + W — > annihilation with probability Pa 



(1.1) 

As such this reaction-diffusion problem has a meaning for any real value < pc < 1, that 
is for any real value q > 2. The Ising case {q = 2) corresponds to a pure annihilation 
problem [p^ = and pa = 1), whereas the q —>■ oo limit corresponds to a pure coagulation 
case {pc = 1 and Pa = 0). It turns out that the later case is much simpler to study than 
any finite q case. In particular, in the q = oo limit, simple random walk arguments are 
sufficient to obtain the distribution of domain sizes 0] and the exponent for persistent 



spins |T^, whereas the computation of the corresponding quantities for any finite q is much 
more involved 0-0]. The relative simplicity of coagulation models {A + A — > A), with 
possibly the back reaction {A — > ^4-1-74) or a random input of A-particles, or localized 
sources of A-particles, is in fact related to the possibility to write closed kinetic equations 
for the "one-empty-interval probabilities", i.e. the probabilities that a given interval contains 



no A-particle |jTT]-0,^ . This approach may be generalized to write closed kinetic equations 



for the probabilities to have many disconnected empty regions ||T^J^, but all these many- 
empty-interval probabilities may in fact be expressed in terms of the one-empty-interval 
probabilities alone f^. This means that all the information on the coagulation model is 
actually contained in these one-empty-interval probabilities. 

The method that has been followed to study the general q case has been to relate 

the zero-temperature Glauber dynamics of the g-state Potts model to another reaction- 
diffusion problem which is a pure coagulation problem (A -|- A — > A) for any q, in contrast 
with the react ion- diffusion model of domain- walls W ( |1.1[ ). A simple way to implement the 
zero-temperature Glauber dynamics consists in updating the spins according to 

(1.2) 
(1.3) 
(1.4) 



S^{t) 


= 5i-l(t 


- dt) 


with probability 


dt 


S^{t) 


= 'S'j+i(t 


- dt) 


with probability 


dt 


Siit) 


= s,it- 


- dt) 


with probability 


(1 - 2dt) 
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Tracing back in time the sequence of spins responsible for the value Si{t) of spin i at time 
t therefore defines a random walk going backwards in time and leading to some spin Sj{0) 
of the initial condition. If we are interested in the values Si{t) and Si{t') of the same spin 
at two different times p| [H, we have to consider the corresponding two random walks 
starting respectively at time t and t' from site i and going backwards in time, and study 
whether they merge at some point before the initial time t = in which case Si{t) = Si{t'), 
or whether they do not meet up to time t, in which case Si{t) = Si{t') with probability ^ 
due to the randomness of the initial condition. This is the starting point for studying the 
probability that a given spin i does not flip up to time t [^]. The same type of reasoning can 
be applied to compare the values Si{t) and Sj{t) of two different spins at time t to lead to 
the distribution of domain sizes 0. In this approach, the Glauber dynamics of the g-state 
Potts model is thus formulated in terms of random walks going backwards in time that 
merge whenever they meet (A + A — > A), and the parameter q only appears through the 
property that two different sites of the random initial condition have the same value with 
probability -. The problem has now therefore a mathematical meaning for any g > 1, in 
contrast with the initial reaction-diffusion model of domain-walls W (|1.1| ) defined for q > 2 
only. This model can moreover be given a physical meaning for any real value g > 1 if one 
considers the Ising case with a random initial condition presenting a non-zero magnetization 
m G [—1,-1-1], in which each spin has initially the value (+) with probability p+ = ^-y^ 



and the value (— ) with probability p_ = |T9|-^,0. Indeed the probability that a (+) 
spin does not flip up to time t or the distribution of (+) domain sizes will be given by the 
corresponding results for the previous g-state Potts model with the correspondence 

This article is devoted to the following generalization of the problem of the exponent of 
persistent spins for the zero-temperature Glauber dynamics of the g-state Potts model. We 
consider a B-particle that diffuses with a diffusion constant Db = XD^ that is different from 
the diffusion constant of domain walls (A G [0, +cx3)), and that disappears whenever it 
meets a domain wall W 

B + W — > with probability 1 

The question is : what is the exponent 6{q, A) that describes the survival probability of the 
B-particle 

PB(t,g,A) « (1.6) 

t — »oo 

as a function of g > 1 and A > ? The exponent is expected to depend explicitly on the 
ratio A = because of the interplay between the diffusion of B and the domain coarsening 
of the kinetic Potts model. From the point of view of reaction-diffusion models, the problem 
considered here is "an impurity problem" , in which there is a small concentration of 

S-particles in an homogeneous background of VT-particles, so that one can neglect reactions 
among impurities and the influence of impurities-background reactions on the background 
properties. In this language, the problem of the fraction of persistent spins in the kinetic 



Potts model can be reformulated as a "static impurity problem" [29 
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What do we know about the exponent 6{q,X)7 It is clear from the definition of the 
model that 6{q, A) is an increasing function of q at fixed A, and an increasing function of A 
at fixed q. For q = I there is no domain walls W so that the exponent vanishes in this limit 
6{q = 1, A) = 0. For a fixed B-particle (A = 0), the survival probability of the B-particle 
reduces to the probability that a given spin is not crossed by any domain wall up to time 
t, which is also the probability that a given spin does not fiip up to time t, and the exact 
expression recently obtained for this exponent reads 



2 2/2 



^(g, A = 0) = -- + ^arccosM (1.7) 



TT 



1] 



V2qJ 



in agreement with previous numerical results |T0|,|23,p|. The exponent 9{q, A) is also known 



in the particular case q = 2 and A = 1 where the i?-particle can be considered as a domain 
wall W, and where the dynamics of domain walls reduces to a pure annihilation model 
{W + W — > 0) [|5]-|2^ 



% = 2,A = l) = i (1.8) 

We have not been able to get an exact expression of the exponent 6{q, A) in the general 
case g > 1 and A > 0, but we have obtained various asymptotic behaviors. The paper is 
organized as follows. In Section H we recall how to derive the value of the exponent 6{q, A) 
in the particularly simple case q = 00 and arbitrary A 



% = cx),A) = - (1.9) 

2 arccos ( y+a j 

In section p| , we extend the approach described in reference by Derrida, Hakim and 
Pasquier to obtain the first correction in e = g — 1 of the exponent for any A 



% = l + 6,A)^e| A +o{e) (1.10) 

^TT - arccos (^YTaJ/ 

As previously explained, this exponent characterizes the decay of the probability that a B- 
particle remains in a (+) domain up to time t for the Glauber dynamics of the Ising model 
starting from a random initial condition characterized by a strong magnetization m = 1 — 2e 
( 11.51) . In section l^, we generalize the approach developed in reference [|] by Derrida, Hakim 
and Pasquier to study the first-order perturbation in A around the result ( |1.7| ) 

e{q,X)=9{q,0) + X4>{q)+o{X) (1.11) 



but the expression obtained for 0(g) is unfortunately quite complicated (see eq ( [4.34|) below). 



To make the reading easier, we gather here the useful notations used in the various parts 
of the paper 



a = a(A)= arccos 



A 



1 + A. 
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/i = fi{q)= — — 
6 = 6{q) = — arccos (— 4/i) 

TT 

6 = 6(q) = — arccos f — ^— , 
TT \V2qJ 



II. DIRECT STUDY OF THE EXPONENT e{q,X) FOR q = oo AND ARBITRARY A 

For q = oo, the dynamics of domain walls (|1 . 1| ) reduces to a pure coagulation model 

W + W — >W with probability 1 

The problem of the B-particle survival is therefore reduced to a three-body problem, since the 
two nearest domain walls enclosing B evolve only by diffusion, and cannot disappear when 
meeting the next nearest domain walls. Let us introduce the joint probability ip{xi,x,X2,'t) 
that the B-particle has not yet disappeared at time t and is at position x, with the next 
domain wall to its left being at position xi < x and the next domain wall to its right being 
at position X2 > x. This joint probability evolves in time according to the diffusion equation 
(where, for simplicity, we have set the diffusion constant of domain walls Dw equal to 1) 

dip d'^ip d'^ip ^ d'^ip „ , , 

^ = TT^ + TT^ + A— ^ for Xi<X<X2 (2.1) 

at oxf 0x2 OX'' 

with the absorbing boundary conditions ■ip{x,x,X2,t) = = ip{xi,x,x,t). The survival 
probability of the B-particle 

/+00 rx r+oo 

dx / dxi / dx2 'ip{xi, X, X2,t) (2.2) 
-oo J—oo Jx 

exhibits the asymptotic behavior (see Appendix A) 

^ / A \ 

-Pb(^; ? = oo. A) oc t 2a(A) where a(A) = arccos (2.3) 

We therefore recover the value of the exponent ( and references therein) 



The angle a(A) decreases from a(A = 0) = | to «(A = oo) = 0, so that 9{q = oo. A) grows 
from 9{q = oo, A = 0) = 1 to 9{q = oo, A = oo) = oo. 
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III. EXPONENT e{q = l + e, A) AT FIRST ORDER IN e FOR ANY A 



In this section, we follow the approach described in reference by Derrida, Hakim 
and Pasquier and only mention the modifications that have to be made for the case we are 
interested in here. 



A. Equivalence with a coagulation model on a large ring 

Finite-size scaling arguments imply that the exponent 6{q, A) also appears in the zero- 
temperature Glauber dynamics of the g-state Potts model defined on a ring of finite length 
L. It describes in this case the decay with the size L of the probability for the B-particle to 
survive indefinitely 

VsiL, q, A) ~ PB(t ~ L\ q, A) ~ L'^^^^'") (3.1) 

As explained in when the spin- values seen by the B-particle are traced back in time, 
one obtains random walks going backwards in time (that we will call A-particles from now 
on), that merge whenever they meet {A + A — > A), and that connect all the spin values 
seen by the B-particle to various ancestors belonging to the random initial configuration. 
The probability for m ancestors of the random initial condition to have the same "color" of 

/ , \ (m— 1) 

the g-state Potts model is simply l-j . As a consequence, the survival probability may 
be expressed as 

rj,(L,q,X)=j2-^PLim,X) (3.2) 

■m=i y 

where ^^(m. A) is the probability of finding m particles on a ring of L sites in the steady- 
state of the following one-species A-particle coagulation problem : there is a moving "source" 
(corresponding to the B-particle) that is always occupied by a A-particle, and the (L-1) other 
sites may be either occupied or empty. During each infinitesimal time-step dt, each A-particle 
hops with probability dt to its right neighbor and with probability dt to its left neighbor, and 
does not move with probability (1 — 2dt), in correspondence with the updating rules of the 
T = Glauber dynamics (p..4[). If two particles occupy the same site, they instantaneously 
coagulate (A + A — > A). Whenever the A-particle being on the source jumps to one of 
its neighbors, a new A-particle is instantaneously produced on the source. The additional 
rules for the dynamics of the source are the following : during each infinitesimal time step 
dt, the source hops with probability Xdt to its right neighbor and with probability Xdt to its 
left neighbor, and does not move with probability (1 — 2Xdt). Whenever the source moves, 
the A-particle that was occupying the position of the source remains on this site, and a new 
A-particle is instantaneously created at the new position of the source. 



To compute the expression (|3.2|) that involves the probabilities ^^(m. A) characterizing 
the steady state of the coagulation model, we introduce generalized "empty-interval prob- 
abilities" |in]-|l8],|l. We first define the conditional probabilities b\f{t) , (l < i < j < L), 
that the segment {s{t) +i, s{t) +j — 1} contains no A-particle, for a given source trajectory 
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{S} = {s{t),t > 0} representing a particular realization of the source random walk. They 
evolve in time according to 



bj^it + dt) = b^^'it) + dt Ssit+<it)Mt) di(^) + btiA^) + bl%,it) + biy_,{t) - Abiyit) 



■{•5} 



{■5} 



l{S} 



+dt 5, 



{t+dt),s{t)+l 



+ dt 6, 



{t+dt),s{t)-l 



and satisfy the boundary conditions 

The average of b\ j {t) over the random walk trajectories of the source denoted by 

{s}. 



i?,,(A,t) = <6,7(t)> 



(3.3) 



(3.4) 



(3.5) 



evolve in time according to 



at 



+A 



-Bj+i,j+i(A, t) + i?i_ij_i(A, t) — 2Bij{\, t) 



(3.6) 



and converge at large time towards stationary probabilities Bij{X) that are solutions of 
Bi+iji\) + 5,_i,,(A) + 5ij+i(A) + 5,,,_i(A) - 45,,,- (A) 



+A 



-Bi+ij+i(A) + 5i_ij_i(A) — 2Bij{\) 



(3.7) 







together with the boundary conditions i?o,j(A) = = i?j^L+i(A) and Bi i{X) = 1 

We may also define the conditional probabilities ^j\/2,---,j2n-i,«2n(^) ? — "^i < "^2 < 
• • ■ < "^2^ < L), that there is no A-particle in any of the segments {s(t) + ii, s(t) + 12 — 
1}, . . . , {s(t) + i2n-iis{t) + — 1} for a given source trajectory {S} = {s(r),r > 0}. 
These conditional probabilities satisfy evolution equations analogous to (|3.3|), with obvious 
boundary conditions for coinciding indices. We have for example 



r 

Kjkli^ + = ^ijkli^) + ^s{t+dt),s{t) 



l{S} 



+dt 5, 



s(t+dt),s(t)+l 



dt S 



s{t+dt),s{t)-l 



^i+ij+l,fc+l,«+l(^) b}j^j^i{t) 

(3.8) 
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The averages < blfj^ ...i^^___^ i^^{t) > over the reahzations of the source trajectories con- 
verge at large time to stationary probabihties -Bn,i2,--,j2n-i,j2n('^) satisfying equations gener- 
ahzing ( p.7|) , as for example 

-Bi+ij,fc,i(A) + -Bi_ij,fe,/(A) + Bij^i^k,ii^) + -Bj,i-i,fc,z(A) 
+Bi 



-85i,,(A) + A Bi+ i,j+i,fc+i,/+i(A) + -Bi-ij-i,fc-i,i_i(A) — 2Bij^k,i{^) 







The generalization of the identity (10) of gives the survival probability of the B-particle 
in terms of the whole hierarchy of the mean empty-interval probabilities -Bji,j2,---,i2„-i,i2n(-^) 
as 

1 



VB{L,q,X) 



qL-i 



1+ Yl {q-iy-'B,,{X)+ Yl (g - l)^-^+'-'Si,,,M(A) + . . . 

^<i<j<L l<i<j<k<l<L 

(3.10) 



The key to solve the coagulation model in the case of a fixed source |^ (corresponding to 
the A = case here) was the possibility to write the mean many-empty-interval probabil- 
ities -Bii,i2, -,j2n-i,i2n('^ = 0) arbitrary order as Pfaffians of the mean one-empty-interval 
probabilities Bi j{X = 0) alone, as for example 

B^,,,k,liO) = 5,,,(0)5,/0) + 5,/0)5,,,(0) - 5,,fe(0)5,/0) (3.11) 

For the case of a moving source, all these Pfaffians relations still hold for a given real- 
ization {S} of the source trajectory. It is for example easy to check that 

^%,iit) = blfit) bifit) + blfit) b^Ht) - blfit) 6jf (t) (3.12) 

since the right hand side and the left hand side evolve in time according to the same equation 
( p.8| ) and satisfy the same boundary conditions. However, these Pfaffians relations that hold 
for a given realization of the source trajectory are no longer valid for the mean probabilities 
Bii,i2,---,i2n-i,i2nWy because Pfaffians involve products and the mean of a product is of course 
not equal to the product of means. In particular, equation ( p.ll|) is no longer valid for A 7^ 

5„-fc/A) ^ fi,,,(A)5fc/A) + 5,/A)%fc(A)-5,,fc(A)%(A) (3.13) 

Equation ( p.lO|) giving the survival probability of the B-particle in terms of the whole hi- 
erarchy of the mean many-empty-interval probabilities -Bji,j2,--,j2n-i.i2n(-^) is therefore much 
more difficult to use for A 7^ 0. 

This is why in the following, we only compute the mean one-empty-interval probabilities 
Bij{X) in the limit of a large system L —>■ 00, and use the following expansion of the survival 
probability (|3.10| ) in e = g — 1 

VB{L,q = l + e,X) = l + e ^ (5,,,+i(A) - l) + 0(6^) (3.14) 

i=i ^ ' 

to obtain the first order in e of the exponent 6'(g = 1 + e. A). 
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B. Mean one-empty-interval probabilities Bij{X) for a large system L ^ oo 

For large L, Bi j{X) becomes a continuous function Px{x = x,y = i) which satisfies the 
continuous version of ( |3.7| ) 

(1 + -^) yg^2 dy'i j + dxdy ^ (3.15) 

in the triangle < x < y < 1 together with the boundary conditions P\{0, y) = = Px{x, 1) 
and f3\{x, x) = 1. To eliminate the non-diagonal term, we perform the change of coordinates 

X = x and y=(i±%^ (3.16) 

V2XTT 

The problem is now reduced to solving the Laplace equation 

^ + ^-0 (317) 

in the triangle of vertices 0{Xo = 0, Fq = 0), A {^Xa = 0,Ya = -^^fj and 
C {^Xc = l,Yc = 755^) 5 with angles OAC = a = arccos(i^) and AGO = CO A = 7 = 
see Figure 1) with the boundary conditions 

(3x{X, Y) = on segments OA and AC and (3x{X, Y) = 1 on segment OC (3.18) 

A convenient way to solve this problem is to use methods of complex analysis and con- 
formal transformations. We introduce the function (3{w) of the complex variable w = u + iv 



TV— a 
2 



P{w) = Q 



1 , /w 
— In 



TC \ w 



(3.19) 



where denotes the imaginary part. This function l3{w) satisfies the Laplace equation on 
the upper half-plane {v > 0} and the following boundary conditions on the real axis {v = 0} 

f3(u e {-00, 0), V = 0^ =0 = p(^ueil,+oo),v = 0^ and G (0, 1), f = 0^ = 1 

(3.20) 

We now consider the conformal transformation 

Z(w) = K / di— — e^(f (3.21) 



with the notations 



'' = ^1 + -) A-^ ^^fP-^-' (3,22) 

2 V TT/ r^(l - v) cos(|) 

This transformation maps the upper half-plane of the complex plane w into the interior 
of the triangle OAC described above in the complex plane Z = X + lY . We have more 
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precisely Zq = Z{w = 0), Zc = Z{w = 1) and Za = Z{w = ±00). The function Px{X,Y) 
of ( |3.17| - |3.18| ) therefore reads 



Px{Z = X + lY) = P[w{Z)] = Q 



(3.23) 



where w{Z) is the inverse mapping of Z{w) (|3.22| ). Unfortunately, the inverse mapping w{Z) 
cannot be explicitly written for arbitrary Z in the triangle. However, approximate explicit 
forms may be written locally. 



C. Use of Bi^l{X) to recover 6{q = 00, A) 



The exponent 6(q = 00, A) ( |2.4| ) may be recovered from the asymptotic behavior of 
since ( p.2|) reduces for g = 00 to 



Vb{L, g = 00, A) = A) 



(3.24) 



and the probability of having only exactly one particle on the ring is simply the probability 
-Bi,l(A) of having no particle except on the source. So finally 



PB(L,g = 00, A) = 5i,l(A) ~ /5A(a; ~ -,?/ ^ ^ 



(3.25) 



We thus have to consider the function (3\{Z) for Z near the vertex A of the triangle. In the 
neighborhood of Z^, w{Z = Z^ + re^^~^~^'^^) reads approximatively for G [0, a] and small 
enough r 



w[Z = ZA + re'^'^^'*''^ 

and we get 
/?,(Z = Z^ + re^(-f+^))^ 
We finally obtain 



) - - 


'a r ' 







(3.26) 



nw{Z = Z^ + re^(~f+'^)) 



1 


■ ar ' 








sm vr— 


TT 


[nK] 


V a 



y ~ 1 - i) ~ P,{Z ^Za + ^e^(-i+<^)) oc L-S (3.28) 
and so recover again the exponent (|2.4|). 



D. Exponent 6'(g = 1 + e, A) at first order in e 

To obtain the first order in e of the exponent 

e{q = l + t,\) = ea{\) + o{e) 



(3.29) 
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we only have to consider the Bi^i^i{X) since a(A) is given by the leading behavior (|3.14|) 



L-l 



^(l-B,.,+i(A))^^ 2a(A)lnL 



i=l 



(3.30) 



We may locally invert Z{w) around Zb = Z{0) and obtain w{Z = re**- 2 '^'^) for (j) G [0, 7] 
and for small enough r 



w{Z = re*(5-'^)) ~ - 
1 — >o 



—I- 



. 



We therefore get 

f3x{Z = re''-^-^^) - ^ 



In 



1 — z/)r 



e 1-" 



1 - 



7r(l — z/) 7r(l — z/) 



(3.31) 



arctand) ^^^^^^ 



or more explicitly in terms of the original coordinates {x, y) 

r, / ^ -, \ 2 / X sin a 
Pa(x -C 1, y <^ 1) ~ arctan I 



TT — a \y — X cos a 



(3.33) 



As explained in this small corner < x < y < 1 where j3x{x,y) is of the form /a(^) 
with the scaling function 



2 / Msma 

jxiu) = arctan 

TT — a \1 — u cos a 



(3.34) 



and the symmetric corner l — y<l — X'^l, are entirely responsible of the leading behavior 
of ( |3.30|) , that is more explicitly 



X:(l-i?v+i(A)),^ 2L fdx 



I - (3x[x,x + 



(3.35) 



=.2L l^da; ^2/j,(l)lnL 
The correction a(A) (|3.30|) is therefore given by 

a(A) = /;(1) 

so that finally 



(vr — a) tan(|^ 



%=1 + 6,A) = 6 I ^^^^^T^ I + 0(e) 



TT — arccos 



l+A 



(3.36) 



(3.37) 



(3.38) 
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IV. EXPONENT 0{q,X) AT FIRST ORDER IN A FOR ANY q 

In this section, we follow the approach of reference |^ by Derrida, Hakim and Pasquier 
and only mention the modifications that have to be made for the present study. 



A. Adaptation of the formalism introduced in [6] 

Let us introduce the survival probability (t, q) up to time t of the B-particle for 
a given trajectory S = {s(r),r > 0} of the B-particle. We may also define the analog 
R^^\t,q) of this survival probability for the semi-infinite-chain geometry [Q. The survival 
probability on the infinite chain reads then as a generalization of Eq (10) of 

p^f{t,q) = R\,^'\t,q)Ri-'Ht,q) (4.1) 

where the notation {—5} denotes the mirror-trajectory {— s(r),r > 0} of the trajectory 
S = {s(r), r > 0}. The generalization of Eqs (29-30-31) of [H is that the survival probability 
Rb (tiyh, q) between times ti and ^2 for a given trajectory S may be expressed as 

R^i\hM,q) = Q^'\h,t,,q) ei^^^>(*-*-^) (4.2) 

where, using the notation = ^ and denoting by di the derivative with respect to the 
first variable of any function of several variables, we have 



n Jti Jti 



n=l 



(4.3) 



and 



with 



Q^''Hh,t2,q) = ^l-fici^}{h,h,q) - {q-l)^-fici^}{t2,h,q) (4.4) 



h,q) = E(-2/i)" f'dr,--- f'dr^ 

n=l -^il 



(4.5) 



So the fundamental object needed to compute R^^\ti,t2,q) is the probability c^'^\ti,T2) 
(with < Ti < T2) that, for a given trajectory S = {s(r),0 < r < T2} of the source, two 
Brownian walkers going backwards in time and starting respectively at s{t2) at time U2 = 
t — T2 and at s{ti) at time Mi = t — ri do not meet up to time t, with the boundary condition 
that they are reflected by the reversed-time source trajectory {<j{u) = s(t—u) , t—T2 < u <t} 
(see Figure 2). 

The method of images gives that c^'^^{ti, T2) may be written as 
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c'f'^^(Ti, T2) = / dx dy 

Ja{t)=s{0) Jx 



(4.6) 

g^^\x,t\ s{T2),t-T2) g^^Hy,t\ s{n),t-Ti) 



in terms of the probability density g^'^^{x,t\ s{T),t — r) that, for a given trajectory S of 
the source, a Brownian walker going backwards in time and starting at s(r) = a{t — r) 
at time {t — r) is at site x at time t, with the boundary condition that it is reflected by 
the reversed-time source trajectory {a{u) = s{t — u)}. More exphcitly, g^^^{x,t\ a{u),u) 
satisfies the diffusion equation 

^^{^>(x, t\ a{u),u) = ^9^^\x, t\ a{u),u) for x > a{t) (4.7) 
together with the initial condition 

g^^^ix, t\ (t{u),u) — > 5[x - aiu)] (4.8) 

and the reflection condition at a; = a{t) expressed by the conservation of probability 
^ r J 15}/ .1 / A ^ n \s\, u\ .1 r ^ n , dg^^^x,t\ a{u),u) 



x=a{t) 

(4.9) 



Unfortunately, we have not been able to write this probability density g^'^^{x,t\ a{u),u) 
in an explicit simple way as a functional of the source trajectory {S}. As a consequence, the 
average over the sources trajectories that is needed to evaluate the exponent 9{q, A) through 



PB{t, q, A) =< p^^\t, q) >=< R\t^\t, q) R'^^^^t, g) > « r^(«'^) (4.10) 



t— *oo 



seems quite difficult to study. 

So from now on, we will restrict ourselves to the evaluation of the B-survival probability 
at the flrst order in the diffusion coefficient A of the B-trajectories. 

B. Perturbation theory in the diffusion constant A 

The probability density g^'^^{a{t) + z,t\ a{u),u) may be seen as the continuous limit of 
a discretized version involving [N — 1) intermediate times tn — u-\- < n < N — 1) 

g^^^a{t) + z,t\ a{u),u) = 

N~l (4.11) 

roc roc roo I— r r i 

lim / dzi I dz2--- dzjq-i TT g^ ^(cr(t„+i) + cr(t„) 

~l~ 1 tn) 

N-^ooJo Jo Jo 

with the conventions to = u, Zq = 0, = t and zj^ = z. In the limit of a vanishing time 
interval At — tn+i — tn — ^ — > 0, we may approximate the source trajectory between 
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times tn and by a line segment of slope a„ = "'('^n+ih-'^iin) _ ^j^g probability density 
g{a(u)=ao+au} 2;o,0) solution of ( |4 ■ 7| - |4 ■ 8| - |4!9| ) for the particular case of a linear trajectory 
cr(u) = ao + au reads 



2V7rt 



4t _^ ea(a;o - ctq) 



4t 



v/ — oo 



(x + Xq - 2o-o - riY 
At 



(4.12) 



So we get 



5f^'^^(0-(t) + Cr(M),M) = lim / rf^i / dZ2--- dZN-l Y\ QaA^n+l, Zn, At) (4.13) 



where 



Qa„{Zn+l, Zn, At) = Q 



{a{u)=a{tn )+anu} , 



aitn+l) + Zn+l,tn + At\ (j{tn) + Zn, tn) (4.14) 



e 4 / 

2yH 



I^^n+l ~ ^n)^ (^n+1 + ^n) 



+a„ I drin e 2 e 4At 



4At +e 



4At 



We then obtain the following expansion up to second order in {s{u)} 

g^^Hait) + z, t\ ait-T),t-T)= go{z, r) + gl^\z, r) + gi^\z, r) + 
with the term of order 



iz,T) 



Z 

e 4r 



vrr 



the term of order 1 in {s{u)} 



(4.15) 



(4.16) 



(4.17) 



,P>(z,r) = -l rdu^^^fe-^ 



^ -^0 ^M(r - u) 



(4.18) 



and the average of the term of order 2 in {(j{u)} 
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( _t\ 

I 4r 



V 



(4.19) 



Note that this last average may be obtained directly from the average < g^^^{a{t) + z^ t\ a{t — 
T),t — t) > which is simple for arbitrary A (See Appendix B). We then obtain the corre- 
sponding expansion 



(4.20) 



c^^Hn, r2) = Coin, n) + c[^\n, n) + ^(n, r2) + ■ ■ ■ 

of c^^^Ti,T2) defined in (|4^ ) 

g^^Ha{t) + z, t\ ait -n),t- n) g^^Ha{t) + z\ t\ a{t - ra), t 



c{'5}(ri,r2) = dz I dz' 



CO roo 



-g^^^ait) + z,t\ a{t-T2),t-T2) g^^\a{t) + z\t\ a{t-n),t-n) 



with the term of order 



4 /rT 

Co(ti, T2) = 1 arctan J — 

TT V '^2 



(4.21) 



(4.22) 



the term of order 1 in {s{u)} 



Cp^(n,T2) = ^ 



7r2 



r2 



du ^^^^ I 

T2 + M y u{ti — u) Jo Ti + f y v{t2 — v) 



•^2 siv) 

dv 



(4.23) 



and finally the average of the term of order 2 in {s{u)} 



<cP(ri,r2) >= 



2A ri du 



"^2 dv 



mm -u, f 



V — u 



{v + uY 



-^0 ^Ju{ti - u) -^0 ^v^r2 - v) 
To compute the first correction in A of the exponent 

e(g,A) =0(g,O) + A0(g) + o(A) 
we have to study the survival of the B-particle between times ti and t2 ( [4.1CI| - ^^ 



(4.24) 



(4.25) 



oc t 



-e{q,\) 



(4.26) 



For 1 < g < 2, the prefactor (5^'^^(ti, ^2, g) Q^"'^^(ti, ^2, g) remains finite in the limit to oo 
(see P for the detailed study of the A = case). Using the expansion 
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(4.27) 



we find that the correction ([4.25| ) is given by the leading behavior of the average 



<Tt^\t^M,q)> -A0(g)ln(-^ 

^2^00 \tl 



(4.28) 



This average at first order in A of T^^^^{ti^t2^q) defined in ( [4. 3D decomposes into two 
contributions 



< Tt^\h,t2, q) >= X(%{h,t2, q) + ri,i(ti, t2, g)) 



(4.29) 



with 



T2(ti,t2,g) = -E(-2/^)"/ c?ri---/ dr^- <d,ci^\n,T2)>Y[dMn,n+i) (4.30) 

n=l -^^l -^^l ^ i=2 



and 



°° ('—2/7')" " /•*2 /•t2 



(4.31) 



SO that the correction (j){q) in (|4.25|) will be obtained through the asymptotic behavior 

r2(ti,t2,g) + Ti,i(ti,t2,g)^^ -0(g) In f^) (4.32) 



The asymptotic behaviors of T2(ti,t2,q) and Ti{ti,t2,q) are studied respectively in Ap- 
pendix C and Appendix D and we only give here the final result for 0(g) obtained in 
Appendix E in terms of the auxiliary variable 



\ 2 /2-g\ 
olg) = — arccos — ^— 
\V2qJ 

which varies in the interval S(q = 1) = 1/2 < 6 < 5{q = oo) = 3/2 



(4.33) 



cos(7r(5) / fT 1 

^ — 1 + 35 + 5^ - + - 

7r2sin(7r(5) V V4; 2 
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4-*jni-*)-(^*)*(^^ 



+- 



1 /-i dz 



vr' Jo V-2 



1(2^'^) ^ ^/ 1 (2 + 



(4.34) 
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where = is the logarithmic derivative of the Gamma function, and where the 

functions H^{z,6) and H^{z,6) are defined as the series 



2 cos(7r(5) 
sin(7r(5) 



n=l 



(2n+l-(5) 



r (2n + 



1-^ 



V Z i ^y^'o I i I V z 

r (2n + 1 - 5) fo r (2n + 1 + 5 



(2n+l+5) 



- rf2n 



1 + 6 



(4.35) 



H_{z,S) = 2 



r 00 
J2 ^{2n+l-5) 

n=l 



T(2n + 2-6 

r (2n + 1 - 5) ' ^ 







{2n+l+<5) 



1^(2^^ + 2 + 5 



r (2n + I + <5) 



00 

/=0 



(4.36) 



The function 0(g) is plotted on Figure 3 in terms of the auxiliary variable 6{q) 
I arccos (^^) in the interval 6{q = 1) = 1/2 < 5 < 6{q = 00) = 3/2. 
It is easy to check that 



and 



0(g = 00) = - 
TT 



0(g = 1 + e) = ( - - 4) e + o(e) 



(4.37) 



(4.38) 



which are consistent with the previous results ( p. 41) and ( 3.38| ) at first order in A 



V. CONCLUSION 



We have studied in this paper a two-species reaction diffusion system in the limit where 
the minority B-species has a very low concentration in comparison with the majority W- 
species, with the following two-particle reactions : two W-particles either coagulate (W + 
W — > W) or annihilate (W^ + iy — > 0) with the respective probabilities pc = (Q'~2)/(g — 1) 
and Pa = l/(g — 1); a B-particle and a W-particle annihilate {W + B — > 0) with probability 
1. We have seen why the exponent 9 {q, X = Db/Dw) describing the asymptotic time decay 
of the minority B-species concentration could be viewed as a generalization of the exponent of 
persistent spins in the zero-temperature Glauber dynamics of the ID Potts model starting 
from a random initial condition. We have extended the methods introduced by Derrida, 
Hakim and Pasquier for the problem of persistent spins |^ to compute the exponent 
6{q, A) in perturbation at first order in (g — 1) for arbitrary A and at first order in A for 
arbitrary q. Let us now briefiy outline the problems that have to be overcome to go beyond 
the first order perturbation theories presented here. 
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In the approach of section |T|, the survival probabihty of the B-particle is given in terms of 
the whole hierarchy of the mean empty-interval probabilities -Bn,i2, --,i2„-i,j2n(-^) of some one- 
species coagulation model with a randomly moving source ( |3.1CI| ). Here we only computed 
the mean one-empty- interval probabilities Bij{X), and we were thus limited to the first order 
in e = g — 1. However, we have seen that for a given realization of the source trajectory, 
many-empty-interval probabilities of arbitrary order could be written as Pfaffians of the 

r c\ _______ 

one-empty-interval probabilities b}j (t) alone ( |3.12| ). One could therefore think of writing 
the survival probability of the B-particle for a given realization of the source trajectory as 
the square root of some determinant, as in formula (11) of reference The remaining 
problem then consists in evaluating the average of the square root of this determinant over 
the realizations of the source trajectories. 

In the approach of section 0, the fundamental object involved in the expression of the 
survival probability of the B-particle for a given realization of the source trajectory ([4.1[ - [4.6|) 
is the probability density g^'^^{x, t\ a{t — T),t — T) that, for a given trajectory S of the source, 
a Brownian walker going backwards in time and starting at a{t ~ r) at time {t — r) is at 
site X at time t, with the boundary condition that it is reflected by the source trajectory 
{a{u)} ([4.7| - |4.9| ). Here we only computed this probability density g^'^^ in perturbation up 
to second order in the source trajectory, and we thus had to restrict ourselves to the first 
order in the diffusion coefficient A. However we have seen that the probability density 
g^'^^(x,t\ o{t — T),t — r) could be written as some functional of the source trajectory {S} 
( [4.13| - |4T5| ) and one could perhaps hope to put this functional in a sufficiently simple form 
to be able to evaluate the exponent 9{q, A). 
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APPENDIX A: ASYMPTOTIC BEHAVIOR OF Psit, q = oo,X) 



To obtain the survival probability Psit^q = oo,X) (^.21), it is convenient to introduce 
the relative coordinates yi and j/2 together with the center of mass g for the diffusion prob- 
lem (13) 

yi=x-xi>Q ; y2=X2-x>Q ; g = 2\ + l ^ ' 

The survival probability then reads 

r+oo r+oo 

PB{t,(X),X)= dyi dy2 f{yi,y2,t) (A2) 
Jo Jo 

where the function / satisfies 

+ -TTT T i for Z/i > and y2 > (A3) 



1 + \ dt dyf dyl 1 + A dyidy2 

with the boundary conditions f{0,y2,t) = = /(?/i,0,t). To eliminate the non-diagonal 
term, we set 

X = y, and y = ^^l±£±i)^ (A4) 

V2XTI ^ 

The problem is now reduced to the diffusion equation 

+ — ^ in the corner Y > —== X > (A5) 



1 + A 9t 9X2 ' gyi V2A + 1 

So finally in polar coordinates, we get 

n+OO 



/'+00 /"Q 

Ps(t, oo,A)=/ pdp\ d(f) g{p,(f),t) (A6) 
Jo Jo 

where g{p,(j),t) is the solution of the diffusion equation 



I dg 1 d f d'ip\ 1 d'^ijj 



p-^ I H — 2^7a corner < < a = arccos ( y+a) {^'^) 



1 + \ dt pdp \ dp ) p 

with the absorbing boundary conditions g{p, 0, t) = = g{p, a, t). The solution g{p, 0, t) for 
a given initial condition 

g{p, 0, t) — . 1 5{p - po) 5(0 - 0o) (A8) 
may be expanded onto an eigenfunction basis as 

g{p, 0, t) = £ r k dk hk,miP, 0) hk,rniPo, 0o) e^^^ + ^)^'^ (A9) 
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with the eigenfunctions 

hk,m{p,(t>) = \ - Jm^{kp) sin (mTT- ) (AlO) 
V q; " \ a J 

where J,y{z) denotes the Bessel function of index p. The survival probabihty of the B-particle 
may finaUy be written as 

PB{t,oc,X)= Y: —sm(m7r^] H r dr T q dq e'i^ + J^^qr) J^J^] 
nTtdd""^ \ OiJJo Jo " "\ViJ 

(All) 

The asymptotic behavior of Bessel function at small argument, J^{z) ~ r{i+u) i^/'^T ^ 
2; — > 0, gives the long-time behavior 



PB{t,q^ 00, A) ^ -sin (tt^) ^ ^ ( /° I oc (A12) 



t-ooTT V ay r(l + ^) V2^(l + A)i, 



i— >oo 



APPENDIX B: DIRECT EVALUATION OF Vx{z,t) =< g^^Ha(t) + z,t\ a(0),0) > 



The mean probability density R\{a, x, t) that the source starting at ctq = at time t — 
is at position a at time t and that the random walker emitted by the source at time t — 
is at position x > a at time t is solution of the diffusion equation 



d_ 
dt 



Rxia,x,t) = lX— + — \Rx{a,x,t) 



ioT X > a 



(Bl) 



together with the initial condition 

and the reflection boundary condition ai x — a 



_ dRx{a,x,t) _ ^^dRxia, x, t) 



dx 



da 







(B2) 



(B3) 



This boundary condition which may be obtained from the continuous limit of the discrete- 
space dynamics, ensures that the partial law for the position of the source is a free Brownian 
motion of coefficient A as it should 



/ dx Rx{o', X, t) — 

J a 



a 

e~4At 



(B4) 
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as can be checked by taking the time-derivative of both sides. 

The mean probabihty density Vx{z,t) =< g^'^^{a{t) + z,t\ a{0),0) > of the relative 
coordinate z = x — a therefore satisfies the diffusion equation 



dVx{z,t) d^Vx{z,t) 



dt 



together with the initial condition 



for 2; > 



Vxiz,t)^5{z) 

and the following boundary condition at 2; = 

dVx{z,t 



dz 



z=0 



So it simply reads 



(B5) 



(B6) 



(B7) 



'Px{z,t) 



4(1 + A)t 
7r(l + \)t 



(B8) 



from which we directly recover equation (|4.19|). 



APPENDIX C: ASYMPTOTIC BEHAVIOR OF T2{ti,t2,q) 
To compute 72(^1,^2,0'), we use the notations 



Co(^i, ^.+1) = go f 1 and < cf^{Ti, T2) >= A 5(2 ( — ) (CI) 



with (|4:22|) 



and ( IC^ ) 



Ji+l ^^2 



4 

gQ{z) = 1 arctanv^ (C2) 

TT 



2 /-i dx dy y-zx 



and rewrite 



%{h,t2,q) = -f:{-2f,r Tdn-.. TdT^-g'J^) n ^^^of^l (C4) 
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/"ill t2 /'in 12 

The leading behavior of order ^ may thus be obtained as in 

g) ^ - In U 5: (-2/x)" / ■ ■ ■ / g'^ (e^i) J] (e^^*) 5{v, +V2 + ---Vn) 

^''l^n=l "^^"^ "^""^ i=2 



Setting 



J-oc 27r 



(C6) 



(C7) 



we obtain 



l2{h,t2,q) ~ - In 



t2\ r+°° dk 



tl 



2ti 



dv e'k^ e g',{en 



n=l 



+ 00 



+ 00 



2^1 + 2/i/o(fc) 



(C8) 



(C9) 



where 

Io{k) -- 
and 

Using (|C3|) , we get 



w 

dw e'kw = _ 2 ^^fct. ^ ^ 2 

TT J-oo 1 + cosh[7Tk) 



+00 



(CIO) 



(Cll) 



h{k) = -Ik- / rfx ^ / dy ^ / dw e^^^min(e-, 1) ,^ ^ " , (C12) 

The integrals over the variables x and y give 

.^fciny r (1 - z/c) r (i + ^A;) _ TT sinh(7rA;) 



1 Q—iklnx «i 

^/a;(l-a;)^o 



TT- 



(C13) 



/^(^l _ ^) r(l-zfc) r(l + ifc) fccosh(7rfc) 
and the integral over the variable w may be rewritten as 

rfu; e'''^min(e"', 1) = -2t / rfMsin(fcM)- — ^ (C14) 
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to obtain 



4 rccj I 

hik) = tanhfvrA;) / dusm(ku)- — (C15) 

TT JO (e" + 1)^ 

We therefore get 

r,(i,i..,) ^ ^^in(^i) ^"^"'f' m (C16) 

TT Vti/J-oo cosn(7rA;) — 4/i 



Vti/io (e" + l)2 7o ^ ^ cosh(7rA;) - 4/i ^ ^ 



The integral over k gives (in the sense of distributions theory) 

f^"^ • /, \ sinhfvr/c) cosh(u6) . ^1 / . x r -, n 

/ aksmiku) r-r~r\ ^ = ^ where o = — arccos(— 4/ij G (U.5, Ij 

Jo cosh(7rA;) — 4yU sinhw vr 

(C18) 

so that 

^/ X 2u, /t2\ /""^ , « coshfw^) ,^ , 

^ Vti/Jo (coshf 



APPENDIX D: ASYMPTOTIC BEHAVIOR OF Ti{ti,t2,q) 
To compute Ti i{ti,t2,q), we introduce the notation ( ^4.23|) 



/ 



n Ta T2 



T2 n n 



I <cf\n,r2)cf\r,,n)> 



(Dl) 



dx 



v/x(l-x) -^0 ^yil-y) 



dy 



^/T2nmm{TiX,Tay) ^ y/TYT^mm{T2X,uy) 

(r2 + Tix) {n + TaX) (n + T2X) (r„ + Ux) 



./r[nmm{T2X,Tay) ^/r^ min^TiX , uy) 



(ri + T2X) {n + TaX) {t2 + Tix) {Ta + Ux) 



and rewrite 



(D2) 



n=2 2 



ar„ (9i(92/ — , , — 



(D3) 
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- E E / dui--- dundid^f (e"l-"^ e-'-'+i, e"^-"'+i) [] d'o (e"-«»+i 



n=2 2 ;=2-^l°*l 



The leading behavior of order ^In ^ j therefore reads 

T,,,it„ t,, ,) ^ - In (I) E 5^ . . . J_ 



+00 



dVr, 



n=2 



i=2 



Setting 



J— 00 ZTT 



and 



(D4) 



(D5) 



J— 00 \ / J— 00 J— 00 2,71 



we obtain 



h{k) 



n—l r 



/-2 



(D6) 



(D7) 



(D8) 



where 



^ .^1/ 27r 



J(A:,p) 



27r 



1 + 2/x7o(A;) 



1 + 2//7o(p) 



(D9) 



J{k,p)= [ "duf "dv [ "d^e^^^ + ^^^^ + ^(^-^^)^e^ + ^5i52/(e«,e^,e^) (DIO) 
J— 00 J— 00 J— 00 

Since 



1 dx 



f 

Jo 



dy 



+ 



^x{l - x) ^y{l-y) 
min(a;, e~^y) min(a;, e'"'~^y) 



6 2 min(e"a;, e'" ^y) 
{l^e^x){l^e"x) 

min(e"x, e'^^y) 



(e2+e 2j;)(e2+e ■^y) {e^ + e 2x){l + e'"x) (1 + e"x)(e2 + e 2^) 
it is convenient to begin with the integration over the variable z 

F (e", e^ k-p)= dz e^(^ ' P> f (e^, e^, e^) 
J —00 ^ ' 



(Dll) 



(D12) 
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and to compute J{k,p) through 



hoo /■+00 



oo J — oo 



J{k,p) = -kp I du I dv e^^^ + 'P"" F (e", k - p) 



Using successively 



dz e 



i{k — p)z 



6 2 minfa, 6e 



ah — 



,i(fc-p)ln(^) 



+ {k-pf 



(D13) 



(D14) 



du 5r 



ae2+6e 2 va6cosh(7rcj 



(for a > and 6 > 0) 



(D15) 



dx 



ik \nx 



TX- 



T{\ + ik 



'0 ^x(l-x) ■ r(l + ^A;) 

and well known properties of the F function, we finally get 

1 



Jik.p) = - 



{k — pY cosh(7rA;) cosh(7rp) 



.cosh(7rA;) 
so that 



sinh(7rp) + sinh(7rA;) ^ 



cosh(7rp) 



— J / dk I dp- 



T{l-ik) T{l-ip) 
cosh(7rA;) cosh(7rp) 



(D16) 



(D17) 



27r^ \tiJ J-00 J -00 cosh(7r/c) — 4yU cosh(7rp) — 4/i 
It is convenient to decouple the integrations over the variable k and p by setting 

1 



J{k,p) (D18) 



\ + {k-pY 



dx e"^ e^('=-P)^ 
J —00 



(D19) 



in (D17) to obtain 



Ti,i(ti,t2,g) ^ -4^1n(^^^ r'^dxe-^ 



\ / J —00 



-00 /■+00 



00 J —00 



-ipx 



dk I dp 



sinh(7rj9) + sinh(7rfc) 



cosh(7rA;) 



kp 



cosh(7rA;) — 4/i cosh(7rp) — 4/i 



cosh(7rp) I r(l — ik) T{1 — ip) 



+ c.c. 



(D20) 



Using again ( P18|) and 



25 



n+OO 



dk sin (/ex) 



k 



1 d 



cosh(7rA;)(cosh(7rA;) — Afi) 4/i dx Jo 



1 d 



dk cos{kx) 



sinh(5x) 



^cosh(7rA;) cosh(7rA;) — 4yU^ 



Ajjidx l2cosh(|) sin(7r5) sinhx^ 



(D21) 



leads to 



1 r°° 
tiJ llJo 



dx e ^ 



cosh(5x) 
cosh^ f 



1 



sin(7r5) Jo 



dx e 2 {26 



\tiJ Jo 



with 



h{x) = — i 
2 



cosh^((5a;) sinh(25x) coshx' 
sinh^ X sinh^ x 

dx e~2 h{x)h{—x) 

T{l-zk 



(D22) 



+CX) t.y:jikx 

dk 

-oo cosh(7r/c) — 4/i r(l — ik) 



d r°° ^ 1 sinh 5(x + z) 
dz- 



^/^^sm{'^■6) dx Jo -\/e^ — 1 sinh(x + z 
For a; > 0, may be easily expanded as 



h{x) 



n oo 

sin(7r5) 



-^(2n+l-(5) 



i-.).r(2n+|-5 



-(2n+l+<5) 



^ r (2n + I + (5 



r (2n + 1 - 5) r (2n + 1 + 5) 

and may be rewritten in terms of hypergeometric functions 6, c, z) ( pO| ) 



(D23) 
(D24) 

(D25) 



h{x) 



sin(7r5) 



r(i-5) 



F ( 1, - - 5, 1 - 5, + ( 1, - - 5, 1 - 5, -e-" 



-'2 ' 



r(i + 5) 



F 1, - + (5, 1 + 5, e-" + F 1, - + (5, 1 + 5, -e 



The analytic continuation of the hypergeometric functions 
sion for h{—x) in the domain x > 



(D26) 

gives the following expan- 



h{-x) 



cos(7r5) 



E 

n=0 



g (2n+i-5). r (2n + 2 - ,5) ^ ^_(2„+i+5). r (2n + 2 + ^) 



r (2n + f - 5 



T{2n + I + S 



(D27) 



26 



APPENDIX E: EXPRESSION OF 



The asymptotic behaviors of 72(^1,^2,9) Ti{ti,t2,q) obtained respectively in Ap- 
pendix B and Appendix C imply that 0(g) defined in ( f4 .251 - ^321) is given for 1 < g < 2 in 
terms of 5 = ^ arccos(— 4/i) = ^ arccos(— 4^^) G (0.5, 1) as 



2 cos(7r(5) 
TT^ sin(7r5) Jo 

cos^(7r5) /""o 



°° , _x ( ^cosh^f^x) sinh(25x) coshx' 
dx e \25- 



TC^ JO 



sinh^ X sinh^ x 

dx e^2 h{x)h{~x) 



(El) 



where the function h{x) is given in ( P26| ) of Appendix C. 

The first integral of the formula ( |E1| ) can be computed in terms of the logarithmic 



derivative \E'(x) = of the Gamma function 



°° , _£ / ^cosh (6x) smh (25a;) cosh X 

dx e ^ { 26 — ^ — 

\ sinh X sinh x 



,1/4 



il-z) 



(5(1 - z) [z^ + z'^ + 2) + (1 + ;2) - z"^) 



(E2) 



3. K 1 
-5+ - + - 

2 2 V4/ 4 



l-.).(i-.)-(i..)*(i.. 



To obtain the analytical continuation of (j){q) to the domain g > 2, we need to redefine 



5 as 



o(gj = — arccos — 4- 



TT 



g- 



^ 2 /2-g\ 
olg) = — arccos — ^— 



(E3) 



and to rewrite ( |E1| ) in a form valid in the whole domain 6{q = 1) = 1/2 < 5 < 6{q = 00) 
3/2 



^, , cos(7r(5) / /1\ 1 



1 /•! rf^ 



i7+(^,^)^-(^,^) + 2;.^'-'^ f^^T^^^ 



r(i-5) 



4z,6)- 5-- 



2^r(i + 5) 



(E4) 



where 



2 cos(7r(5) 
sin(7r5) 



1-5 



n=l 



r (2n + I + 5 

r(2n + l-(5) ^0'' r(2n+l + 5) 



_ ^ ^(2n+l+<5) 



(E5) 



27 



■n=l 



(2n+l-(5) 



r{2n + 2-6 
r{2n+l-S) 

tA^H) £(i±l) 



n=0 



(2n+l+<5) 



r(2n + 2 + S) 
r (2n + I + ^) 



1=0 
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